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Abstract: We derive some exact results concerning the anomalons U(1 )a symmetry in 
the chirally symmetric phase of QCD at high temperature. We discnss the importance of 
topology and finite-volnme effects on the U(1 )a symmetry violation characterized by the 
difference of chiral snsceptibilities. In particnlar, we present a reliable method to measnre 
the anomaly strength in lattice simnlations with fixed topology. We also derive new spectral 
snm rnles and a novel Banks-Casher-type relation. Throngh onr spectral analysis we ar¬ 
rive at a simple alternative proof of the Aoki-Fnkaya-Tanignchi “theorem” on the effective 
restoration of the U(1 )a symmetry at high temperatnre. 
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1 Introduction 

The physics of the U(l) axial symmetry in quantum chromodynamics (QCD) has been a 
subject of intensive research over many years. While the QCD Lagrangian is invariant under 
X \ J ( Nf)L at the classical level, the U(l)yi symmetry is broken by quantum effects 
[1] (see (2.1)). In fact, experimentally observed hadron spectra in the vacuum do not fully 
respect the U(Aj)/j x \J{Nf)L symmetry: the r]' meson concerning the U(l)yi symmetry is 
much heavier than the other pseudoscalar mesons associated with chiral symmetry breaking. 
This U(1 )a problem was resolved [2, 3] through the discovery of nonperturbative topological 
excitations in the Euclidean spacetime, called instantons. 

Recently, much attention has been focused on another “U(l)y 4 puzzle” in QCD at 
T > Tc, where Tc is the (pseudo)critical temperature of the chiral transition. Namely, 
even though the U(1 )a anomaly relation (2.1) itself does not receive any modification at 
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finite temperature [4]/ the U(1 )a symmetry could be effectively restored at the level of 
mesonic two-point (or higher-point) correlation functions [5]. This problem is of particular 
interest due to the role played by the axial anomaly in determining the order of chiral phase 
transition [6]. It was claimed by Cohen [7] that massless two-flavor QCD at T > should 
be effectively symmetric under U(l)y 4 in the sense that two-point correlators in the vr, a, 5, 
and r] channels become all degenerate. This work was soon followed by counterarguments 
[8, 9]. Later it was recognized by Laine and Vepsalainen [10] that the U(l)yi violation in 
the flavor-singlet (axial) vector channel can be shown rigorously at least for high enough 
temperatures without any assumptions. This conclusion qualitatively agrees with [11] that 
shows U(l)yi violation at high temperatures by calculating a nonzero splitting of scalar 
and pseudoscalar screening masses using a semiclassical dilute instanton gas picture [12]. 
More recently, Aoki et al. [13] claimed to have shown the effective restoration of the U(1 )a 
symmetry rigorously under certain assumptions, although the validity of their assumptions 
appears to be nontrivial and subtle (see Sec. 5). 

Alongside these theoretical studies, the U(l)y 4 problem at finite temperature has also 
been studied intensively in first-principles lattice QCD simulations [14, 15], but a consensus 
is not reached yet: effective restoration of the U(1 )a symmetry was reported in simulations 
with overlap fermions [16] and domain-wall fermions [17-20], whereas a violation of the 
U(l)yi symmetry was reported in simulations using staggered fermions [21-24] and domain- 
wall fermions [25].^ We warn that some of the simulations [17, 21-25] were performed for 
2 -|- 1 flavors; the effect of a heavier strange quark on the possible U(l)y 4 violation in the 
light-quark sector is not completely clear yet. 

In this paper, we do not try to solve this U(1 )a puzzle at finite temperature. Rather, 
we derive some rigorous results involving the U(1 )a symmetry in high-temperature QCD, 
assuming that the U(1 )a symmetry is violated. (Since the U(1 )a violation must be present 
at least for high enough temperatures [10, 11], we consider its presence at all T > Tc 
to be quite plausible.) We first derive general expressions for chiral susceptibilities and 
the topological susceptibility at T > Tc using the method of [15, 26]. They are used to 
highlight that, while the majority of the U(1 )a violation at small volume comes from exact 
(topological) zero modes, the dominant contribution at large volume comes from nonzero 
modes. We estimate finite-volume effects suffered by lattice QCD simulations with fixed 
topology, and propose a way to measure the U(l) a- violating effects reliably in a finite 
volume. Furthermore, we rigorously derive new sum rules and a new Banks-Casher-type 
relation for the Dirac eigenvalue spectra at T > Tc- This relation provides a link between 
the connected two-point correlation function of Dirac eigenvalues and the U(1 )a anomaly. 
As a by-product of our spectral analysis, we find a remarkably simple proof of the Aoki- 
Fukaya-Taniguchi “theorem” on the effective restoration of the U(1 )a symmetry at high 
temperature, under the same assumptions as in [13]. It should be emphasized that all 

^We noticed that the tensor decomposition of the anomalous correlation function at finite temperature in 
[4] is incomplete, and that many more terms need to be included; see Appendix A. This does not, however, 
affect the conclusion of [4] that the anomaly relation (2.1) remains unchanged at finite temperature. 

^In [24] the overlap Dirac operator was used to probe the Dirac spectra while configurations were 
generated with improved staggered fermions. 
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our results are based on a systematic analysis of QCD. We expect that testing our exact 
relations and proposal in future lattice simulations should be a useful step towards the 
resolution of the U(l)yi puzzle at finite temperature. 

The paper is organized as follows. In Sec. 2, we review the argument of [10] for U(1 )a 
violation at high temperature. In Sec. 3, we discuss the importance of topology and finite- 
volume effects on the breaking of the U(1 )a symmetry, as well as its implication for lattice 
QCD simulations. In Sec. 4, we derive new spectral sum rules and a Banks-Casher-type 
relation for Dirac spectra concerning the U(l)^ anomaly. In Sec. 5, we comment on the 
Aoki-Fukaya-Taniguchi “theorem” in [13]. Section 6 contains our conclusions. 

In Appendix A, we point out and correct the deficiency in the tensor decomposition of 
the anomalous correlation function in QCD at finite temperature studied by Itoyama and 
Mueller [4]. In Appendix B, we discuss the microscopic scaling that is different from (4.12) 
in the main text. In Appendix C we derive (4.19) in the main text. Throughout this paper, 
we will work on QCD with quarks in the fundamental representation of the gauge group 
SU(3). 

2 U(l)^ anomaly at high temperature 

In this section we review an argument for the U(l)yi anomaly in massless QCD at high 
temperature, given by Laine and Vepsalainen [10].^ Their argument is based on the anomaly 
relation 

6 ^ 1 -'“ = ( 2 , 1 ) 

for the axial current and the Debye screening of the gauge fields at high 

temperature. Here Njis the number of flavors, g is the QCD coupling constant and is 
the gluon field strength with a being the color index. 

Let us work in Euclidean spacetime with the imaginary time r = it G [0, /?] and a 
spatial box of size Li x L 2 x L 3 . We assume that both quarks and gauge fields obey 
periodic boundary conditions in spatial directions, whilst quarks (gauge fields) obey the 
anti-periodic (periodic) boundary condition in the temporal direction. As our interest is in 
the screening of gauge fields, we shall consider a spatial correlation function. Without loss 
of generality we choose a spatial separation in the 0:3 direction. Then, for the axial “charge” 
(integrated over the volume transverse to the X 3 axis) 

Qtixs) ^ J , (2.2) 

one can show that {y^)^ = 0 in the thermodynamic limit if the axial symmetry is 

effectively restored, and 7 ^ 0 if it is broken. Here (...) is a statistical aver¬ 

age, and the coordinates x^ and ^3 are entirely arbitrary. The fact that {Qi{x^)Q'i{y'i)) = 0 
when U(l)/i is unbroken may be shown in two steps as follows. Suppose we add a constant 

®We thank M. Laine for explaining the logic of [10] in detail. 


- 3 - 




external field a that couples to the axial charge Q^. Then the action acquires an additional 
contribution 6S = a dxs Qfix^). This leads to the susceptibility 


A 1 *9 


V 4 da^ ^ Qtixs) j dys (2.3) 

= J dx3 (Qs (x3)Q3 (0)) , (2.4) 


where V 4 = /IL 1 L 2 T 3 is the total volume of spacetime. As U(1 )a is conserved by assump¬ 
tion, must be finite in the thermodynamic limit; hence < 00 as L 3 —)> 00 . This 
completes the first step of the proof. 

In the second step, we use the fact that Q^i^s) is independent of X3 . Actually, this is 
a direct consequence of the (assumed) conservation of the axial current: 


^<3^(a;3) = J drcfx^dsj^ 

= - / dr d^x± {dojQ + V± ■ jf) = 0 , 


(2.5) 

( 2 . 6 ) 


where the last step follows from the boundary conditions for fields. Note that (2.6) holds 
rigorously in a finite volume. Then it is clear that we must have 


^ {Qi{x3)QiiO)) = ^^<33 (a^3)Q^(0)^ = 0 . 


(2.7) 


Combining (2.7) with the finiteness of x^ in limit L 3 —>• 00 , we conclude that (Q'^(a:: 3 )Q^( 0 )) 
must vanish in the thermodynamic limit for an arbitrary X 3 . This completes the second 
step of the proof. 

Next, recall that the anomaly relation (2.1) can be rewritten as a total derivative, 

d,j^ = d,K„ K, = . (2.8) 

In terms of this K^, the so-called Chern-Simons charge reads as 

QP(a;3) = J (ird^x^Ks. (2.9) 

Because the gauge fields are Debye screened at high temperature, the correlation of Qg ® 
should decay exponentially: 


(QP(^3)QP(2/3))~e-l"3-s/3l/« 


( 2 . 10 ) 


at 1 x 3 — 2/31 where ^ is the correlation length (or the inverse Debye mass); ~ gT 
at leading order in g. Although the correlator (2.10) may appear to vanish trivially due to 
the non-gauge-invariance of Q 3 ®, this is not necessarily true. As Qg® is gauge invariant up 
to surface terms, it becomes gauge invariant and its correlator becomes nonzero once we fix 
boundary conditions for the gauge fields. 
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Since and Qg ® have the same quantum numbers, Qa ^ contributes to the correlator 
of Q^. Hence we expect 

{Qiix3)Qtiy3)) OC (QP(x3)QP(y3)) ~ e"l"^3-i/3|/« ^ ( 2 .11) 

This shows that {Q^{x^)Q^{y^)) cannot be a constant as a function of X 3 , indicating that 
the U(1 )a symmetry is certainly violated at the level of correlation functions of quark 
bilinears in QCD at sufficiently high temperature. 

It would then be quite natural to expect, by continuity, that the U(1 )a symmetry be 
violated at any T > Tc- In the following, we shall assume U(l)yi violation for T > Tc in 
the thermodynamic limit and pursue its consequences in detail. 

3 Topology and finite-volume effect 

In this section we write down the most general QCD partition function for T > Tc in terms 
of quark masses and derive general expressions for the chiral susceptibilities and topological 
susceptibility based on the method of [15, 26]. Our arguments here are based on symmetries 
of QCD and a systematic expansion in terms of a small parameter m/T <C 1, and are fully 
under theoretical control. We then elucidate the contributions of topology and finite-volume 
effects to the violation of the U(l)yi symmetry (characterized by the difference of two-point 
functions x-w ~ XS to be defined below), and discuss possible implications for lattice QCD 
simulations. For definiteness, we will concentrate on two-flavor QCD below. 

3.1 Partition function and topological susceptibility 

We consider the partition function of two-flavor QCD at finite temperature as a function 
of quark masses Since there are no massless modes at T > Tc in the chiral limit, ^ the 

free energy density should be analytic in quark mass.^ 

To write down the general form of the free energy, we consider a generic quark mass 
matrix M (=2x2 matrix in the flavor space) and let M transform under the symmetry, 
Q = SU(2)/j X SU(2)j;, X U(1)a, so that the quark mass term in the QCD Lagrangian 

Cmass = IplMlpR + h.C. (3.1) 

is invariant under Q. Here V'it.L ^re the right- and left-handed quarks, which transform 
under Q as i/jr —)> &'‘^^Vr'4)r and ipi —)> respectively, where 6a is a U(1 )a 

rotation angle and Vr^l £ SUr^l( 2). It then follows that M should transform as M —)■ 
e- 2 i(?A VlMVI- 

Noting that the free energy density at T > Tc is invariant under the restored SU(2)/j x 
811(2)/, chiral symmetry but not under the U(l)y 4 symmetry, the partition function of QCD 

“^In the imaginary-time formalism, contributions of massless quarks is infrared (IR) finite because the 
lowest Matsubara frequency for fermions ~ ttT acts as an effective IR cutoff. 

®It should be stressed that analyticity of the partition function breaks down if the zero-temperature part 
is thrown away. For example, in a non-interacting theory the free energy of quarks after subtraction of the 
zero-temperature part includes a term ~ m^log(m/7rr) [27], which is not analytic in m. 


- 5 - 



in a spatial volume V 3 can be expanded in terms of a small parameter mu^/T <C 1 as 
[15, 26] 


Z(T, Fg, M) = exp -^/(T, ^ 3 , M) 
f{T, Vs, M) = fo- f 2 tr MtM - /^(det M + det ikft) + ©(M^), 


(3.2) 

(3.3) 


where /o, /2 and are functions of T and V 3 . We assume that this expansion has a nonzero 
radius of convergence. The term oc Ja represents the effect of axial anomaly: for a U(1 )a 
rotation ip —)• this term transforms as detM —?• det M, so it breaks U(1 )a 

down to Z 4 . The absence of 0{M) terms is consistent with the vanishing chiral condensate 
in the chiral limit for T > Tc- In the following we will disregard the 0{M^) terms in the 
free energy as they are suppressed by additional powers of mu,d/T <C 1. Since the partition 
function (3.2) is obtained with a systematic expansion, this will be called the “effective 
theory” in this paper (although there is no dynamical field in it). 

We now turn to the study of topological sectors. As is well known, the 9 angle can be 
incorporated into the partition function via M —>• [28], where Nf = 2 is of our 

interest here. Then the partition function in a sector of given topological charge 

32^/ (3.4) 

is obtained, from (3.3), as 

Zq{T, Vs, M) = e-'«^ Z{T, Vs, 

= p-V4[fo-f2iml+ml)] I -iQ6> ^ 214 / 4 cos 6> 

J 2tt 

= e-'"d/o-/2(m2+™2)] j^^2V4fAmumd), 

where V 4 = Vg/T is the spacetime volume, Iq is the modified Bessel function of Q-th order, 
and M = d\a,g{mu, rnd) was substituted. Intriguingly, the probability distribution of Q is 
proportional to Iq in one-flavor QCD, too [28].® 

The Taylor expansion of (3.7) in powers of quark masses starts with (y4,fA'mufnd)^^\ 
which is the contribution of exact zero modes. Hence the topological sectors with Q 7 ^ 0 
will all drop out in the chiral limit if V 4 is finite. By contrast, topological fluctuations will 
not be suppressed at all even near the chiral limit if V 4 is sufficiently large. This subtle 
balance between topology and volume has an important practical consequence for lattice 
simulations, as we will discuss shortly. 

An important quantity that characterizes topological fluctuations is the mean square 
of the topological charge at 0 = 0 , 

Q=—oo 

®An analogous toy model was also studied in [29, Appendix A]. 


(3.5) 

(3.6) 

(3.7) 
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where (3.7) was used. The topological susceptibility is then given [15, 26] by 


Xtop — 


Vi 


2fAmumd. 


(3.9) 


Alternatively, one can reach (3.9) by considering that, with the replacement M ^ M 
in (3.3), the 0-dependence of the free energy reads 


/(^) = / - 2 fAmumd cos 6 + 


(3.10) 


where / = /o — f 2 {^u + ^d) is the term independent of 6. The topological susceptibility is 
then 

= 2fAmumd, (3.11) 

9=0 

which is the same as (3.9). Note that our result, obtained assuming Ja 7 ^ 0, does not agree 
with the result by Aoki et al. [13] that {Q'^)/V^ —)• 0 as V 3 —?■ 00 for a small but nonzero m. 
The topological susceptibility here should be contrasted with that of the QCD vacuum, 

(3.12) 

with S being the magnitude of the chiral condensate [28]. This difference can be under¬ 
stood in the following way. In the presence of chiral symmetry breaking, the quark mass 
dependence of the free energy can be expanded in terms of the quark mass as 


Atop — 


d^f{9) 




f = fo- S[tr(M[/^) + tr(Mtt/)] + C>(m2), (3.13) 


where U denotes the SU(2) Nambu-Goldstone field associated with spontaneous chiral sym¬ 
metry breaking. The topological susceptibility is dominated by the second term in (3.13) 
and is given by (3.12), with the higher order contributions being O(M^). When chiral sym¬ 
metry is restored (S = 0 ), on the other hand, the leading contribution to the topological 
susceptibility is O(M^) and is given by (3.9). The behavior xtop oc ijiumd is also found in 
the 2SC phase of dense QCD for the same reason [30]. 

From the partition function (3.3), higher moments of Q can also be derived [26] as 


{Q^) 

= A, 

(3.14a) 

{Q^) 

= A{l + 3A), 

(3.14b) 

iQ^) 

= . 4(1 + 15 ^ + 15 ^ 2 )^ 

(3.14c) 

{Q^) 

= A{1 + 63.4 + 210.42 + 105.4^), 

(3.14d) 


with A = 2VifA'mumd , while all odd moments vanish. 


3.2 Chiral susceptibilities and U(l)^ anomaly 

We now turn to two-point correlation functions of quark bilinears (also called chiral suscep¬ 
tibilities). The utility of chiral susceptibilities in two-flavor QCD as a convenient probe for 
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the U(l)/i anomaly has been advocated long time ago [5, 7], and nowadays they are mea¬ 
sured in lattice simulations with dynamical quarks [16, 17, 25]. It is therefore of primary 
interest to relate these susceptibilities to the coefficients /o, /2 and of the effective theory 
(3.3) [15]. In this paper we will work with the following definitions for chiral susceptibilities: 


Xa = 
Xtt = 
X5 = 

Xri — 


d'^x {ipi'y5T3'ipix)'ilji'y5T3ip{0)) - 


d"^x (V’r3V'(x)V’r3'0(O)) - 


d"^x (V’i75V’(a^)V'n5V'(0)) - (V’^TsV’)^ 


(3.15a) 

(3.15b) 

(3.15c) 

(3.15d) 


where '0V’ = uu + dd and = ui'j^u -|- di'j^d in our notation, and T 3 is the third Pauli 

matrix. Note that these definitions do not necessarily coincide with those in the literature. 
Inserting M = diag(mu,md) into /(T, V 3 ,M), we obtain 


X<7 = 


X5 = 


32 q 2 q 2 

+ 7W^+2- 


V 4 \dm^ dm^ drriudmd 


+ - 2 - 


52 


1 


V 4 \dml dm?, dmudmd 


log Z = 4 / 2 + 4/^ 
log Z = 4 / 2 -4/^ 


(3.16) 

(3.17) 


up to 0(m2) corrections. Similarly, inserting M = diag(mu -|- ib, md — ib) into f{T, V 3 , M) 
we find 


1 

1 d'^Z 


fl dZ 

T 

4^4 

z 952 

b=0 

[z db 

b=oJ 


4/2 + 4/a . 


(3.18) 


In the same way one can also show Xr] = 4/2 — 4/a. The equality Xa = Xtd as well as 
XS = Xrj, is a direct consequence of the restored SU(2)/j x SU(2)i chiral symmetry. 

The axial anomaly manifests itself in the difference [15] 


X- - X 5 = 8 /a + 0 (m 2 ) , 


(3.19) 


where the correction is displayed for completeness. Practically, the correlators of vr and 6 
are most convenient in lattice simulations, because they have no disconnected components 
(for degenerate masses). Note that, if one wishes, it is entirely straightforward to extend the 
calculation for (3.19) to higher orders by including U(l)A-violating terms such as (detM)2 
and (tr MMf )(det M) in the free energy (3.3). However this is not expected to bring about 
quantitative differences when M/T <C 1. 

In [11], the two-point correlators ( 7 r(x) 7 r( 0 )) and ((l(x)<5(0)) (rather than their spatial 
integrals, x-w and Xs) were calculated directly in high-temperature QCD by using the ’t 
Hooft vertex of instantons. They found that the difference of the two correlators decreases 
at high temperature but does not vanish exactly, in agreement with the general argument 
presented in Sec. 2 . 
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It was emphasized in [8, 9] that the dominant contribution to Xn — XS comes from 
exact zero modes in the Q = ±1 sector. A more recent paper [13] argues to the contrary 
that contributions of exact zero modes is suppressed in the thermodynamic limit. In what 
follows we aim to clarify this issue. 

Let us first decompose the anomalous contribution (3.19) into contributions from each 
topological sector. We assume 0 = 0 in the following. Since the second terms in (3.15b) 
and (3.15c) vanish for degenerate masses, it follows that 


lim (xtt -Xs) = / d^x 

J 

1 

oo 

- E 

Q=—oc 


('tpij5T3'lp(x)i’i75T3'lp(0)) - ('lpT3lp(x)tpT3l/^(0)) 


/ 1 d'^Z 


1 d'^Z 

[z 962 

6=0 

Z 9c2 


c=0 


Z 


(3.20) 

(3.21) 


where it is tacitly assumed in (3.20) that the first term is evaluated for M = diag(m + 
ib,m — ib) and the second term for M = diag(m + c,m — c). In (3.21) we defined the 
contribution Pq from the sector of topological charge Q as 




V4 


( 1 d^Zg 


1 d^ZQ 

\Zq 962 

6=0 

Zq 9c2 


4/2 + 4/a 


= 8/a 


lQ{2VifAm^) 
lQ{2VifAm^) ■ 


c=0 


4/2 - 4/a 


I'Q{2V^fAm^ 

lQ{2VifArrP 


(3.22) 

(3.23) 

(3.24) 


Using the identity Iq{x) = ^Iq{x) + Iq+i{x) or I'q{x) = -^Iq{x) + Iq-i{x) depending 
on the sign of Q, one may cast Pq into a suggestive form 



fQ+i(2V4/Am^ 


for Q > 0 , 
for Q < 0 . 


(3.25) 


The first terms in (3.25) are the contributions from exact zero modes. This can be easily 
seen by plugging Zq oc (m^ + 6^) I *31 and Zq oc (m^ — c^)l*3l into the first and the second 
terms in (3.22), respectively. Therefore the U(l)A-violating contribution (3.21) may be split 
into the zero-mode fraction' and the nonzero-mode fraction as 


lim {xn - X5) = 8fA{Sz + Snz) , 


(3.26) 


^It is intriguing that (3.27) below has exactly the same form as the fraction of zero modes for the chiral 
condensate in one-flavor QCD [28, eq. (7.3)], under the identification 2VifA'm7 -e-)- UEm. 
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X 

Figure 1. Relative contributions of zero and nonzero Dirac eigenmodes to Xir ~ Xs ^ ^ function 
of a; = 2 VnfA'm^- 


where 

^ (2 E [lo{2V^fAm^) + 11(2^4 W)] , (3.27) 

Snz = l-S,. (3.28) 

In addition, the contribution of the Q = ±1 sectors to Sz is defined as 

The quantities Sz, Snz and S±i are plotted in Figure 1 as functions of x = 2 V 4 ^fAm^- 

We observe that, in a small volume or near the chiral limit (x <C 1), Xn — X5 is 
dominated by the contribution of exact zero modes in the Q = ±1 sector, as argued in 
[8, 9]. By contrast, if we take the thermodynamic limit (x 1), the contribution of 
nonzero modes dominates, and the exact zero modes are completely irrelevant. This can 
be understood from (3.8): since (Q^) ~ V^fArn^, one naturally expects (|(5|) = 0(\/l4), 
implying that the first term in (3.25) is suppressed in a large volume.® On the other hand, 
the second term in (3.25) tends to 8/a, which is the same value as in the full theory (3.19). 
This means that the anomaly (/a ^ 0) in the thermodynamic limit must be attributed to 
nonzero Dirac eigenmodes. The Q = ±1 sector does not play a distinguished role. Indeed, 
one can show for x 1 that ZqjZ obeys a Gaussian distribution (see also [28]), according 
to which ZqjZ ~ l/Y^V4/Am,2 for |<5| < \/V 4 Ja^ and is suppressed otherwise. Therefore, 
if the volume is sufficiently large with a fixed nonzero mass, all contributions to Xtt ~ XS 
from the sectors with |Q| < -s/VaJa^ are equally important, in contradistinction to the 
finite-volume regime (x < 1) where only the Q = ±1 sectors contribute to Xn — XS ■ 

To avoid confusion, we stress that the total amount of Xtt~X5 is equal to 8/a irrespective 
of the value of x; the order-of-limit issue does not arise, of course, because there is no long- 
range-order in QCD above Tc . The reason the exchange of dominance occurs between zero 

®In this inspection, the positivity of the path-integral measure plays an essential role. We note that the 
suppression of exact zero modes does not hold in general for negative or complex masses [28, 31, 32]. 
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Figure 2. The magnitude of (xtt — X5)|q^o normalized by (xn — X5)|fun ns a function of x = 
2T4 /aw^- At large volume (x » 1), /i(x)//o(x) ~ 1 — 


and nonzero modes as we vary the volume is that a long-range correlation is induced once 
the global topological charge is fixed [28]. 


3.3 Implications for lattice QCD simulations 


We now discuss implications of the above results for lattice QCD simulations. So far the 
U(l)/i anomaly at high temperature has been thoroughly investigated on the lattice (as 
reviewed in Sec. 1), but despite efforts, a definitive conclusion on the (non-)restoration of 
the U(l)/i symmetry is not reached yet. This is not surprising, considering that the physics 
of U(l)^ anomaly is highly sensitive to the explicit breaking of chiral symmetry by lattice 
discretization; even domain-wall fermions have serious problems, as pointed out in [20]. In 
this regard, the most reliable simulations are those in [16] employing dynamical overlap 
fermions. They reported restoration of the U(l)y 4 symmetry based on simulations with 
a fixed global topological charge (Q = 0). They also evaluated possible finite-size effects 
associated with the topology fixing, by using the formalism developed in [29, 33]. Here we 
wish to revisit this issue based on our effective-theory framework. 

It follows from (3.21) that in the topologically trivial sector (Q = 0) we have 


Xtt 


- X<5 = 8/a 


Ii{2V^fAm^) 
WVifAm^) ■ 


(3.30) 


The ratio of (3.30) to (xw — X( 5 )|fuu = 8 /a is plotted in Figure 2. It shows that the ratio 
tends to 0 for small x and obscures the nonzero value in the full theory. This signals a strong 
finite-volume effect at small x. It seems necessary to ensure at least x = 2 V 4 /Am^ > 1 in 
order to observe a nonzero value of Xtt — X& clearly. 

Our result so far is rigorous, as long as /a 7 ^ 0 and the 0{M^) correction to (3.3) 
can be neglected. At sufficiently high temperature T ^ Tc we may resort to the dilute 
instanton gas approximation [ 12 ], which yields 


/a ~ ~ ^ J.-23/3 ('g 3 ^^ 


for iVc = 3 and Nf = 2 . Since it decays so rapidly, it would be a challenging task to achieve 
a sufficiently large volume that satisfies 2 V4fATn? > 1 while keeping m small. On the other 
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hand, near , the asymptotic formula (3.31) breaks down and we do not exactly know how 
small Ja is. 

One way to extract /a from a topology-fixed simulation is as follows: if the simulation 
volume V 4 is not large and hence < 1 , then we have, to a good approximation. 


(Xtt 


X<5) 


<3=0 




(3.32) 


from (3.30), where we used Ii{x)/ Iq{x) ~ x/ 2 for x < 1. Therefore in principle one can 
extract /a by fitting the lattice data to the formula (3.32). This is a proposal for future 
lattice simulations with overlap fermions. 


4 Dirac eigenvalue spectra and U(1 )a anomaly 

In this section we derive new spectral sum rules and a novel Banks-Casher-type relation 
which link the Dirac spectrum to the violation of the U(1)a symmetry in high-temperature 
QCD. We then discuss possible forms of the spectral functions. Throughout this section we 
will focus on two-flavor QCD. We denote the purely imaginary eigenvalues of the Euclidean 
Dirac operator D = 7/^(5^ -|- igA^) by {i\n}n with G IR and define the spectral density 
for a fixed gauge field in a finite volume as 

(4.1) 

n 

Chiral symmetry of the Dirac operator, { 11 , 75 } = 0, ensures p^{\) = 


4.1 Spectral sum rules I: macroscopic limit 

The partition function of two-flavor QCD in the sector of topological charge Q is given by 


^Q = 


V YM 

+ mu )]^(1A£ + md) ) ■ 

k i ' Q 


(4.2) 


where the average is taken with respect to the pure Yang-Mills action, and \Q\ exact zero 
modes are implicitly included in the product. Equating the microscopic partition function 
(4.2) to that of the effective theory (3.7) and taking their derivatives with respect to 
and/or nid, one finds various nontrivial formulas for correlation functions of the Dirac 
eigenvalues. The simplest one is given by 


d 

dmu 


log^Q 



2V4/2m« -h 2VifAmd 


lQ{2VifAmumd) 

lQ{‘2yifAmumd) 


+ 0{m^) , 


(4.3) 

(4.4) 

(4.5) 
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with Iq{x) = dlQ{x)/dx. Now (. . .)q stands for the average with full Nf = 2 QCD measure. 
Ill the thermodynamic limit, the number of Dirac eigenvalues scales linearly with V 4 , so let 
us define the one-point function (or the macroscopic spectral density) by 

in terms of which (4.5) reads as 

[ dX = 2f2mu + 2fAmd + 0{m^). (4.7) 

y_oo iX + mu 

Note that i?i(A) depends on and implicitly through the QCD measure used for 
averaging. i?i(A) is expected to have no dependence on Q since topology is irrelevant once 
the thermodynamic limit is taken. Strictly speaking, we have to specify a UV cutoff scheme 
in order to make (4.7) fully meaningful. Equation (4.7) will be used later in Sec. 5. 

There exists another relation that directly relates Ri{X) to the U(1 )a anomaly [15, 34]. 
Setting mu = m^ = m and evaluating the chiral susceptibilities (3.15) in the basis of Dirac 
eigenstates, one can straightforwardly obtain 

dA }? — m? 

By subtraction we arrive at [15, 34] 

roo TUp' 

Xn-Xs = 8 dA — 2 --- ^Ri{X). (4.9) 

Jq [X^ + m^y 

Combining this formula with (3.19), we find 

foo 2 

/ dA RyX) = fA + Oim^). (4.10) 

Jo (A^ m^y 

It is clear from this relation that small Dirac eigenvalues are necessary for Ja to be nonzero.® 
Indeed /a = 0 follows immediately if i?i(A) has a spectral gap near zero at T > Tc. 
Equations (4.7), (4.9) and (4.10) highlight essential properties of i?i(A). 

Unfortunately the form of Ri{X) in the near-zero region cannot be deduced uniquely 
from (4.7) and (4.10) alone. Actually there are infinitely many functions that satisfy (4.10); 
e.g., 

f ^2n+l-k \k 

Ri{X) ~ — 277 - and /^^^^(A), (4.11) 

(A^ -|- m^y 

where n > 1 and 0 < k < 2n + 1 are arbitrary integers. Recently, the Dirac spectrum 
in two-flavor QCD at T > Tc has been studied intensively in lattice QCD simulations 
[16, 17, 24, 25] (see also [14, 15, 35-39] for early works). The three possibilities Ri{X) ~ m, 
A, and m?‘5{X) were examined in detail in [17, 25], whereas the Breit-Wigner form Ri{X) ~ 


®Large eigenvalues with density i?i(A) ~ A® only affect the 0{mJ) part; Ja receives no contribution. 









PqA/{\^ + A^) was nicely fitted to the lattice data in [24, 40] (but see [19, 20, 41, 42] for 
detailed investigations of lattice artifacts stemming from partial reweighting). The 5 form 
is motivated by the dilute instanton gas picture [12] in QCD at T —)• oo, but the exact 
5 form is unlikely to emerge at T > Tc due to the overlap of neighboring instantons and 
anti-instantons. For the moment, contrasting results from different simulations do not allow 
us to draw a dehnitive conclusion on the form of i?i(A). 


4.2 Spectral sum rules II: microscopic limit 

There is yet another way to take the thermodynamic limit in (4.5): if we let all of A, rriu and 
scale as 'i-j'/VAfA in the V 4 —)> oo limit, the dependence on the topology does persist. 
To see this, let us dehne a rescaled dimensionless spectral density 


PQ{C;^^u,Pd) = ,hm 


V 4 —^■oo A 


c 


^JWATa 


. Mu 

'V^WIa' 


md= 


— 


(4.12) 




which is analogous to the microscopic spectral density in the e-regime [28, 43] but note that 
the relevant scale of eigenvalues here is Ij^/V^ rather than 1/14.^^ Then (4.5) becomes 



- Pq{C, P-u, fJ-d) 

< + Pu 


f2 ^qiPupd) 

— Pu + Pd-f—f T ■ 
IA lQ[PuPd) 


(4.13) 


Notice that all 0{m^) corrections in (4.5) drop out in this limit; hence (4.13) is exact. This, 
of course, comes with a caveat that such a limit is meaningful only when the assumption 
/a 7 ^ 0 is correct. (See Appendix B for another microscopic scaling.) It is straightforward 
to extend the rescaling (4.12) to higher-order spectral correlation functions. 

We conjecture that spectral fluctuations on the scale Xjy/Vi^ should be universal, i.e., 
determined solely by global symmetries and independent of the detailed form of QCD 
interactions in the ultraviolet. Although such a new “microscopic limit” prompts us to 
construct a random matrix theory that describes the Dirac spectrum in this regime, we 
have not been successful yet. The difficulty in hnding a proper random matrix theory may 
have something to do with the fact that no global symmetry is spontaneously broken at 
T > Tc, unlike in the QCD vacuum [28, 43] and high-density QCD [44, 46, 47]. 

We can derive inhnitely many spectral sum rules rigorously, by expanding the following 
expression in powers of quark masses. 


(rUnTUrf) QZq 

lim [{mumd)~^ZQ] 
m^O 


^. WVifAmumd) 

{VifAmumd)^ 



(4.14) 


^°Intriguingly, a similar unusual scaling ~ l/\/T4 also appears in color-superconducting phases of QCD 
at high density [44, 45], in the superfluid phase of two-color QCD [46, 47] and in an exotic phase proposed 
by Stern [48-50]. 


14 - 

















where Q > 0 is assumed and the product runs only over A„ > 0. The average above 
is taken for the massless two-flavor QCD measure. Note that both sides of (4.14) are 
normalized to unity in the chiral limit. Two examples of the sum rules read 



) = ^ 4/2 and 


E' 


_ mA? 
K/q 1+q 


(4.15) 


above denotes a sum over > 0. Note that the first sum rule receives UV-divergent 
contributions from large eigenvalues with density oc A^, implying that /2 is a regularization- 
scheme-dependent quantity. By contrast, the second sum rule is dominated by contributions 
from 0{\/ eigenvalues. UV eigenvalues only give 0(14) correction to the RHS, which 
is negligibly small as compared to the 0(14^) term in the V 4 —)• 00 limit. Indeed one can 
show that /a is free from UV divergences (see Appendix B of [49]). The suppression of the 
second sum rule for large Q is due to the repulsion of Dirac eigenvalues from the origin by 
Q exact zero modes. 

In terms of the “microscopic” spectral density, the sum rules (4.15) read 



pq{C,o,o) 

e 



pq(C;0,o) 


A 

2/a’ 


4(1 -|- Q) 


(4.16a) 

(4.16b) 


To obtain the universal function Pq{Q] 0, 0) analytically is an important open problem that 
deserves further investigation. 


4.3 New Banks—Casher-type relation 

If the axial anomaly is present at high temperature, it will be manifested not only in the 
Dirac eigenvalue density but also in the n-point spectral correlation functions for n> 2. To 
examine this possibility, let us introduce the connected two-point correlation function (see 
e.g., [51])ii 


i?c(A,A')= lim ^ 

V4—^00 1/4 


(/(A)p (A')) - (/(A)) (p^(A')) 


(4.17) 


i?c depends on ruu^d implicitly through the averaging weight. Note that Rq satishes the 
constraint 


dA i?c(A,A') = / dA'R c(A,A') = 0. 


(4.18) 


If eigenvalues are entirely uncorrelated, they obey the Poisson statistics. In this case, the 
two-point function is related to the one-point function (cf. [52, (5.4)] and [53, (3.33)]) as 

{p^Wp^{^'))po ^ ■*“ (P^W) ~ ] v ) ( P ^( AO ) > ( 4 . 19 ) 

the rest of this subsection we will ignore topological zero modes altogether. As argued in Sec. 3.2, 
this is justifiable in the macroscopic limit with a positive path-integral measure. 
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where N denotes the number of chiral pairs of Dirac eigenvalues (hence the to¬ 

tal number of eigenvalues is 2N). The <5-functions in hrst term represent a trivial self¬ 
correlation. In Appendix C we outline the derivation of (4.19) for completeness. Then we 
obtain for the uncorrelated Dirac spectra 

i?g°(A, A') = {h(A - A') + h(A + A')} Ri{X) - ^Ri{X)Ri{X '). (4.20) 

As N grows linearly with V 4 , V^/N has a well-dehned thermodynamic limit. One can check 
that (4.20) satishes (4.18). 

Nontrivial two-level correlations among eigenvalues can be characterized by the devi¬ 
ation of Rq from the Poisson case. Let us dehne the two-point cluster function T 2 {X, X') 
by 


T 2 {X, A') = i?P°(A, A') - Rc{X, A'). (4.21) 

We observe that r 2 (A, A') = T 2 (A', A) and 


T2{X, a') = T2{-X, a') = T2iX, -A') = T2i-X, -A') (4.22) 

owing to chiral symmetry. For the Poisson distribution, T 2 vanishes identically by dehnition. 

To see how Rq and T 2 are related to the axial anomaly, we note that (3.16) and (3.17) 
together with Xa = Xtt imply 

X7T-XS = 4xdisc = 8 fA -h O(m^), (4.23) 


where Xdisc is the disconnected scalar susceptibility dehned by 
1 (92 


Xdisc = ,lim 


logZ 


Vi-s-oo V4 dmudrud 

-((y ——y ^ 

/a \ \ i\- -U m. i 


= lim — 

V4—i-oo 1/4 


dA 


1 


iXi -\-mu^ iXj -I- rud 

j 


y —^^ 

^ iXk -\-mr / \ ^ 7 4- 


iXi + nid 


dA'- 


Rc{X,X') 


{iX + mu){iX' + md) ' 

In the limit rriu = rud = m,, we substitute (4.21) into (4.24), obtaining 


/ 


00 

Xdisc — / dA 
1—00 

fOO 


Ri{X) ^ Ri{X) 


{iX -|- m)2 A2 -|- m2 
T 2 (A,A') 


V 4 

N 


'“dA^ 

lA m 


f 


/ dA / dA'^ ^ 

/_oo J-00 [tX + m){iX'+ m) 

o ^^2 coo coo ^ 

( A2 -K m2)2 (iA -h m) (iA' -F m) 


T 2 (A,A') 


= 4/a- 


fdA fdA V' 

y_oo J-00 (^A-km)(^A'-Km) 


(4.24) 


(4.25) 

+ C>(m2) (4.26) 

(4.27) 
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where (4.7) and (4.10) have been used. Recalling (4.23) we obtain 



T2{X,X') 

{iX + m){iX' + m) 


2/a + 0 { m ^). 


(4.28) 


Finally, taking the chiral limit in (4.28), we arrive at a new Banks-Casher-type relation 

72(0,0) = 4/a (4.29) 

for massless two-flavor QCD at T > Tc. To the best of our knowledge, (4.29) is a new 
result. It reveals that the U(1)a anomaly is encoded, not only in the spectral density as in 
(4.10), but also in the nontrivial two-level correlations among near-zero eigenvalues. If the 
near-zero Dirac eigenvalues are entirely uncorrelated, then T 2 vanishes and leads to /a = 0, 
suggesting effective restoration of the U(1)a symmetry. Whether an analogue of (4.29) can 
be derived in Nf > 2 QCD at T > Tc is an interesting open problem. 

Some remarks are in order. In taking the chiral limit we have replaced l/[{iX+m)(iX' + 
m)] with 7r^(5(A)(I(A'). Strictly speaking, in doing so we have tacitly assumed that the 
typical scale over which T 2 (A, A') varies is much larger than m. Whether this is true or 
not in actual QCD is a dynamical problem and must be checked separately.We also 
remark that (4.29) cannot be extended to T < Tc because of the infrared-singular behavior 

T 2 {X, a') ~ log |A - A'l for |A - A'| « A [51]. 

One might suspect that the correlation in the Dirac spectra revealed by (4.29) is at 
variance with the quasi-instanton picture proposed in [26] where a Poisson distribution of 
topological objects (i.e., dressed instantons called quasi-instantons) was argued at all T > 
Tc In the limit T —?■ 00 , where the interaction is weak, the quasi-instanton gas is expected 
to reduce to the conventional dilute bare instanton gas [12]. Let us try to explain how they 
can be consistent with each other. The point is that the Poisson distribution of topological 
zero modes (quasi-instantons) does not necessarily mean the Poisson distribution of Dirac 
eigenvalues. 

In the quasi-instanton picture, independently distributed topological charges are ex¬ 
pected to generate small Dirac eigenvalues that can be described, to a good approximation, 
by a spectral density with a (5-peak at the origin. Let us discuss how to deal with this case 
explicitly within the present spectral analysis. The spectral density over a gauge field 
now assumes a form 


/(A) = c^,5(A) + p^(A), (4.30) 

where p^(A) is the density of eigenvalues away from zero and > 0 is an integer which 
is equal to the total number of topological objects, N = N-^- + N-X^ If we assume that 

^^This smoothness condition is necessary to derive the original Banks-Casher relation, too [28]. 
^^Topological objects similar to our quasi-instantons have been advocated for color-superconducting 
phases of QCD at high density [30, 54[. While quasi-instantons in hot QCD do not interact with each 
other [26[, those in dense QCD weakly interact via exchange of (pseudo) Nambu-Goldstone modes [30, 54[. 

^^Precisely speaking, the ca modes consist of [Qj = \N+ — N- \ exact zero modes and Min(Af+, N-) chiral 
pairs of near-zero modes. 
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the density of nonzero eigenvalues p^{X) is so small that the anomalous contribution /a in 
(4.10) solely originates from the 6 peak at the origin, then it follows that 


(c^) = 2VifAmumd. 


(4.31) 


(In deriving this we have used dx 5{x) = 1/2.) It is pleasing to see that (4.31) agrees 
with the quasi-instanton density derived in [26, Sec. V]. 

Next, by plugging (4.30) into (4.17) we find 


i?c(A,A') 


Vi +-RcO,A') 

+ ym [{c^p-HX')) - (c-A) (pA(A'))] 

+ ySm [{p-'PW) - (c^> (,i^(A)>] , 


(4.32) 


where Rq is the two-point connected function for eigenvalues away from zero. The last 
two lines will vanish if there is no correlation between zero and nonzero modes, which we 
assume. Then we substitute (4.32) into (4.24) to obtain 


1 


Xdisc — 


14 m„mrf L 






dA 


dA' 


Rc{\\') 


{iX + mu){iX' + rrid) 


(4.33) 


As we have been assuming that the density of nonzero modes is sufficiently low, it follows 
that the second term can be neglected in the chiral limit compared to the first term. Then 
(4.31) and Xdisc = 2 /a (cf. (4.23)) imply 

{{c^f) - {c^f' = (c^) = 2VifAmumd. (4.34) 


This coincidence between the average and the variance of ca indicates that is Poisson 
distributed. This is indeed what the quasi-instanton picture in [26] suggests. 

We mention that the Poisson statistics of topological objects was indeed observed in 
recent lattice data at T = 1.5Tc [24]. However, in the real world, quasi-instantons will not 
be strictly noninteracting (due to the 0{m^) term in the free energy) and the (5-peak of the 
spectral density may not be sufficiently narrow to rigorously justify the above treatment. 
Also, the correlations between zero modes and nonzero modes will not be negligible in 
general. With these caveats in mind, we still believe that the quasi-instanton picture in [26] 
and the exact Banks-Casher-type relation (4.29) can be a useful starting point for a fuller 
analytical and numerical investigation of the Dirac spectrum in QCD at high temperature 
in future. 


5 Comment on the Aoki—Fukaya—Taniguchi theorem 

Contrary to our assumption that /a 7 ^ 0 for T > Tc, Aoki at al. [13] claim that, under 
certain assumptions, the violation of the U(1)a symmetry is invisible in correlation functions 
of scalar and pseudoscalar quark bilinears for T > Tc in two-flavor QCD. (This claim does 
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not generalize to the vector-axial-vector sector, as we discussed in Sec. 2.) There are two 

n c c} im ni-i rtm c in rVioiv QTiQlTrcio*^^ 


key assumptions in their analysis: 


1. The Dirac spectral density can be expanded in Taylor series 



(5.1) 


n=0 


near the origin, with a radius of convergence that does not vanish in the chiral limit. 
In particular, there is no 6 {X) term. The notation {Pn)m makes it clear that these 
coefficients are dependent on the quark mass m. 

2 . The expectation value {0{A))m of any m-independent observable 0{A) is an analytic 
function of m? at T > Tc. (Here 0{A) must be a functional of the gauge field only. 
It does not include fermionic observables, such as the chiral condensate. The quarks 
must be integrated out before this assumption is applied.) 

It should be noted that none of the examples in (4.11) satisfies the first assumption. 

Precisely speaking. Ref. [13] assumes that the spectral density for a given gauge field 
Afj, can be expanded in Taylor series, while the above assumption 1 is only concerned with 
the spectral density averaged over all gauge fields.^® 

While their original proof [13] is rather involved, we now show that a much simpler 
proof of /y 4 = 0 for Nf = 2 based on our analysis in the former sections is possible. Namely, 
one can easily prove the following theorem: 

Theorem. Under the two assumptions above, /a = 0. 

Proof. From (5.1) and the Banks-Casher relation [55], we have 


(X \\Ta. {p^)m ■ 


(5.2) 


For T > Tc where (V’V’) = 0) it must be that 


m—>-0 


lim (po = 0 . 


(5.3) 


Since {pQ)m is an analytic function of m? according to the second assumption, (5.3) means 
that {pQ)m = O(m^); in particular 


'o /m ~ 


lim (Po )m = 0{ml). 

TTlu^V 


(5.4) 


On the other hand, it follows from (4.7) that 



(5.5) 


^®Aoki et al. used overlap fermions on the lattice to regularize UV divergences. This is not crucial in the 


following discussion, however. 

AfV-\T- , 


‘We thank Sinya Aoki for clarifying this point to us. 
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where a UV cutoff A was inserted. Note that, to derive this expression, we have only used (i) 
analyticity of the free energy and (ii) irrelevance of exact zero modes in the thermodynamic 
limit (as explained in Sec. 3). If the limit mu —>• 0 is taken with m^ hxed, the RHS of (5.5) 
converges to 2 fA'md + 0 {m^) 

Next, we introduce an arbitrary scale e > 0 which is smaller than the radius of conver¬ 
gence of (5.1) in the chiral limit. Then we split the LHS of (5.5) as 


dA 


2 m,, 


Ri(A) = j^dX y^^Ri(A) + f^dX 


A^ -I- ml 


A^ -I- ml 


sfdA 

Jo 


A^ -I- ml 


(5.6) 


A^ -I- ml 




. 71=0 


n/m , 

n\ 


+ 


2 m„ 


dXRi{X), (5.7) 


where (5.1) was used. The second term in (5.7) is obviously C)(m„), whereas the first term 
is more nontrivial. To check its behavior near the chiral limit, we use 


i: 

L 


dA 


A 


A^ -|- ml 


= - log 


m„ 


+ 0 {ml) , 


'dA 


vr, 


\2 I 2 

A^ -|- mf. 2 


/'dA 

Jo '' ' "-u 

re yn 




t), 


-.n—1 


dA 


A^ -|- m: 


-|- 0 {ml) for n > 4. 


/o ^1" '"'u n — I 
As the leading term in the limit m„ —>■ 0 comes from (5.8a), we deduce 


(5.8a) 

(5.8b) 

(5.8c) 

(5.8d) 


"dA 


10 A2 -h ml 


. 71=0 


A^ 


n! 


_ r 2 mu 

Jo A2 + m2 


{Po)m + 0{mu logmu) 


= [vr 0 {mu)] {Po)m + 0{mu log mu) ■ 
Plugging this into (5.7), we observe that 

2 m,, 


dA 


A^ -|- m; 


:Rl{X) < T^{Po)m + 0{mu logm„). 


Recalling (5.4), it is now clear that 


(5.9) 

(5.10) 

(5.11) 


lim 

m„— aO 


dA 


2 m,, 


A^ -|- m; 


-R^{X) = 0{ml). 


(5.12) 


This is to be compared with (5.5), which tells that the leading term in the limit m„ —>■ 0 is 
2fAmd- Thus = 0 is concluded. This completes the proof. 


^^The 0{m\) contribution comes from higher-order U(l)A-violating terms (e.g., tr(MM^) detM -|- c.c.) 
in the free energy. 
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This short proof is made possible by treating mu and as two independent variables, 
unlike the original one [13], where only the case of degenerate masses was considered. Of 
course, whether the two assumptions are correct or not in QCD is highly nontrivial. If /a 
is non-vanishing in QCD, which has been assumed in the former sections, then one has to 
relax at least one of the two conditions above. Considering that recent lattice simulations 
[17, 24, 25] have demonstrated a singular peak structure in Ri{X) at small A, it seems 
natural to abandon the Taylor expansion (5.1). This issue deserves further investigation. 

6 Conclusions 

In this paper, we derived some rigorous results on the violation of the U(l)yi symmetry in 
two-flavor QCD at T > Tc , which is characterized by the difference of chiral susceptibilities, 
Xtz — X5 (see (3.19)) and is parametrized by /a (see (3.3) for the definition). We clarified how 
the different topological sectors conspire to violate the U(1 )a symmetry and how it varies 
with the spatial volume of the system. We demonstrated that any moment of the topological 
charge at T > Tc can be obtained, once just a single parameter /a is fixed. We also derived 
new spectral sum rules and a Banks-Casher-type relation that relate the anomaly strength 
Ja to statistical correlations in Dirac spectra. As a by-product of the sum rules, we found 
a simple proof of the Aoki-Fukaya-Taniguchi “theorem” on the effective restoration of the 
U(l)/i symmetry [13]. Since nontrivial assumptions are required to prove this theorem, 
we cannot conclude U(l)yi restoration in QCD yet. However, our simplified proof would 
hopefully serve to understand the importance of these assumptions more clearly. 

All of our new exact relations can, in principle, be tested on the lattice. In particular, 
the relation (3.32) can be used to extract the value of /a at T > Tc even in a small volume 
with fixed topology {Q = 0). Finally, we note that determination of /a should also be 
important from a phenomenological point of view, as it is related, through (3.9), to the 
temperature-dependent mass of the QCD axion—an input for the evolution of the axion 
density, which might account for the dark matter density of the universe; see, e.g., [56, 57] 
for recent works. 
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A Tensor decomposition of anomalous correlators at finite temperature 

Consider the anomalous three-point function in QCD at finite temperature T in momentum 
space, 

T^^>^{q,p,T) = Id^xd^y (j"(x)j^(y)j^^(0)), (A.l) 
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where = ■07^V' is the vector current and = '07^7^';/) is the axial current. The color 
and flavor degrees of freedom are suppressed for simplicity. We take the rest frame of the 
medium as r]^ = (1,0). 

Itoyama and Mueller wrote down 30 tensor invariants composed of q^, p‘', r]'^, or 
for [4]. However, their expression is not complete; one can write down 30 more tensor 
invariants. We And that the most general decomposition for T'^pp is given by 


T^PP = + A 2 Pre^^PP + A^rj^E^^PP 

+A^q<^q^ppE^^PP + A,p^q^ppE^^PP + A^r^^ q^ppE^f^PP 
+AjqPqaPpE°‘^'^P + AspPqaPpE^^f^^'P + Agr]PqaPpE°‘f^'^P 
+Aioq‘^qaP^E'^^PP + Aup^qaqpE'^^PP + Ai 2 p'^qar]pE°'^PP 
+AnqPqarii3E°'f^'"P + AupPqaP/sE'^'^^'P + Ai^pPqa’qpE^'l^'^P 
+AiQq^qaPiiE'^^PP + Anp"^ qaqpE’^^PP + Aisp'"qar]pE°‘^PP 
+AigqPpaqi3E''^‘'P + A20pPpa'ni3E°‘^'"P + 

+ [^ 22 ?“^/ + A23p'^qP + A 2 iiq‘"qP + ^25g'^P^ + A26p''pP 

+A27V‘"p^ + A2 %q''r]P + A29p'"r]P + Hsor/'^r?^ + A^ig^^P] qaPi3rj^E°‘^'^P 
+A32qPqaPp^°‘'^‘'P + AsspPqaPi3E°'^''P + A^^g^q^PfiE^^^^P 
+A35qPqagpE'^^''P + A^QpPqaP^E°‘^''P + 

+A^sqPparifiE^^‘'P + A'igpPpaqpE^'^’^P + AiogPpaqpE^^''P 

+ [A4iqPqP + A42pPqP + A^^gPqP + A44qPpP + A4^pPpP 

+A4QgPpP + A47qPgP + A4spPpP + A4ggPpP + ^ 505 ^^] gaP/3??7e“^'^°' 

+ [A^iq'^qP + A^ 2 P''qP + A^sp'^qP + + A^^p^pP 

+H56r/>^ + A57q'^pP + A^sp'^pP + A^gp'^pP + Aeog^^P] qaPpP'yE^'^'^P (A.2) 

with coefficients Ai = Ai{q^, p^, k^, q ■ p, p ■ p, T). The terms H3i^...^eo are new, which 
have been missed in [ 4 ] but are generally allowed by symmetries. It turns out that H34, 
A37, A40, H43, H46, H48, H49, H53, H57, H58, and A^g are not independent of the others 
appearing above and can be omitted without loss of generality.^® In particular, we need 
not consider tensor invariants whose axial-vector index p is carried by pP. Therefore, there 
are 49 independent invariants in total. 


For example, there exists an identity for the T 46 term, 


p u olQ'so p olQop I / \ p aScTu, . / \ p aBcrp , p a aB'ya 

qaPllP VjV E =-qaPpp e ^ + {p ■ pjqap ppe ^ + (q ■ p)Pi 3 p'^qa£ + qappp p^p £ 


(A.3) 


with which the j 446 term can be expressed in terms of the Tg, T 14 , A 20 , and j 427 terms. 
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B Another microscopic scaling 


In (4.12), all dimensionful quantities are rescaled by On the other hand, it is also 

allowed (from a mathematical point of view) to use -^ 214/2 for rescaling. Defining 


w) - ^ 

we obtain from (4.5) a modified spectral relation 


r,, I ,, , , fA 

/ dC - Pq[C,Hu,M = IJ-U +-rfJ-d -ry-^ 

V-oo ^C + ^^u f2 




' Ua, 


(B.l) 


(B.2) 


We suspect, however, that this scheme is precarious because /2 is dominated by UV- 
divergent contributions from the perturbative Dirac spectra i?i(A) ~ which generally 
depends on the regularization scheme. In contrast, /a is free from UV divergences (see 
Appendix B in [49]). This leads us to consider the microscopic scaling by /a as the most 
natural one. 


C Derivation of (4.19) 


The relation (4.19) for uncorrelated Dirac spectra can be easily shown as follows. With 2N 
Dirac eigenvalues we have, from (4.1), 

/ N N \ 

(/(A)p^(A')) = ( {5(A - Afc) + d(A + Afc)} {5(A' - A,) + ,5(A' + A,)} \ (C.l) 

\fc=i £=i / 

= “ ^k) + <^(A^ + Afc)}^ 

+ /{<^(A — Afc) + 6{X + Afc)} {<5(A^ — A^) + 6{X' + A£)}\ (C.2) 

fc=l \ / 

= |5(A-A') + <i(A + A')}(/(A)) 

+ ~ ^k) + <^(A + Afc)^ /{<5(A^ — Xi) + 6{X' + A^)} \ . (C.3) 

k=l \l^k / 

In the last step we factorized the average, which is justified by the absence of correlations 
among different eigenvalues. Then, using the trivial identity 


/^{5(A'-A,)+J(A' + A,)}\ = ^^(p^(A')) , 


(C.4) 


we arrive at the desired formula 

“ ^0 + <5(-^ + -^0} {p^W) “I— {p^W) {p^{^)) ■ (C-5) 

One can see by integrating over A and A' that both sides are normalized to 4A^ correctly. 
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